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Abstract 

Understanding brain connectivity has become one of the most im¬ 
portant issues in neuroscience. But connectivity data can reflect either 
the functional relationships of the brain activities or the anatomical 
properties between brain areas. Although one should expect a clear 
relationship between both representations it is not straightforward. 
Here we present a formalism that allows for the comparison of struc¬ 
tural (DTI) and functional (fMRI) networks by embedding both in a 
common metric space. In this metric space one can then find for which 
regions the two networks are significantly different. Our methodology 
can be used not only to compare multimodal networks but also to ex¬ 
tract statistically significant aggregated networks of a set of subjects. 
Actually, we use this procedure to aggregate a set of functional (fMRI) 
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networks from different subjects in an aggregated network that is com¬ 
pared with the anatomical (DTI) connectivity. The comparison of the 
aggregated network reveals some features that are not observed when 
the comparison is done with the classical averaged network. 


1 Introduction 

In the last decade, the use of advanced tools deriving from neuroimaging 
and complex networks theory have significantly improved our understanding 
of brain functioning [23], Notably, connectivity-based methods have had a 
prominent role in characterising normal brain organisation as well as alter¬ 
ations due to various brain disorders mu Most of the recent works aim 
to quantify the role of connectivity in the communication abilities of neural 
systems. However, the very same notion of connectivity is controversial since 
data used in brain connectivity studies can reflect functional neural activi¬ 
ties (electrical, magnetic or hemodynamic/metabolic) or anatomical proper¬ 
ties BEJ- Neuroanatomical connectivity is meant as the description of the 
physical connections (axonal projections) between two brain sites pE], whereas 
functional connectivity is defined as the estimated temporal correlation be¬ 
tween spatially distant neurophysiological activities such as electroencephalo- 
graphic (EEG), magnetoencephalographic (MEG) functional magnetic reso¬ 
nance imaging (fMRI) or positron emission tomography (PET) recordings [Tj. 

In recent years, the concept of “brain networks” is becoming fundamental 
in neuroscience BUM. Within this framework, nodes stand for different 
brain regions (e.g. parcelated areas or recording sites) and links indicate the 
presence of an anatomical path between those regions, or a functional depen¬ 
dence between their activities. In the last years, this representation of the 
brain has allowed to visualize and describe its non-trivial topological prop¬ 
erties in a compact and objective way. Nowadays, the use of network-based 
analysis in neuroscience has become essential to quantify brain dysfunctions 
in terms of aberrant reconfiguration of functional brain networks [2j 3j [5j • 

Experimental evidence has revealed, for instance, alterations in functional 
and anatomical brain networks in normal cognitive processes, across develop¬ 
ment, and in a wide range of neurological diseases (see PUB and references 
therein). Despite its evident interplay, comparison of anatomo-functional 
brain networks is not straightforward M. Theoretical studies provide 
support for the idea that structural networks determine some aspects of 
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functional networks [B], but it is less clear how the anatomical connectiv¬ 
ity supports or facilitates the emergence of functional networks. Although 
nodes with similar connection patterns tend to exhibit similar functionality, 
the functionality of an individual neural node is strongly determined by the 
pattern of its interconnections with the rest of the network [7]. 

Correspondences between functional and structural networks remains thus 
an active research area mum- A better understanding of how anatomical 
scaffolds support functional communication of brain activities is necessary to 
better understand normal neural processes, as well as to improve identifica¬ 
tion and prediction of alterations in brain diseases. 

In this paper we address this relationship between anatomical and func¬ 
tional connectivity. In previous studies, the correspondence of these net¬ 
works has been often assessed by the difference in an Euclidean space of 
vectors containing connectivity measures such as the clustering coefficient, 
shortest path length, degree distribution, etc. Here, we propose a radically 
different framework for studying brain connectivity differences. Instead of 
extracting a vector of features for each network (anatomical or functional), 
we jointly embed all of them in a common metric space that allow straight¬ 
forward comparisons. Before embed functional and anatomical networks into 
the common metric space, we aggregate group of subjects (e.g. functional 
networks) according to [TB] to obtain a group representation network. There¬ 
fore the method employed in this work allows to preserve in the aggregation 
connected components and to identify among different subjects, a common 
underlying network structure. Our approach may provide a useful insight for 
the analysis of multiple networks obtained from multiple brain modalities or 
groups (healthy volunteers versus patients, for instance). 

2 Methods and Materials 
2.1 fMRI and DTI data 

In this study we considered anatomical and functional brain connectivities 
(extracted from diffusion-weighted DW-MRI and fMRI data, respectively) 
defined on the same brain regions. Brain images were partitioned into the 90 
anatomical regions (N = 90 nodes of the networks) of the Tzourio-Mazoyer 
brain atlas da using the automated anatomical labeling method. 

The anatomical connectivity network is based on the connectivity matrix 
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obtained by Diffusion Magnetic Resonance Imaging (DW-MRI) data from 20 
healthy participants, as described in m The elements of this matrix repre¬ 
sent the probabilities of connection between the 90 brain regions of interest. 
These probabilities are proportional to the density of axonal fibers between 
different areas, so each element of the matrix represents an approximation of 
the connection strength between the corresponding pair of brain regions. 

The functional brain connectivity was extracted from BOLD fMRI resting 
state recordings obtained as described in [13]. All acquired brain volumes 
were corrected for motion and differences in slice acquisition times using 
the SPM10 software package. All fMRI data sets (segments of 5 minutes 
recorded from healthy subjects) were co-registered to the anatomical data 
set and normalized to the standard MNI (Montreal Neurological Institute) 
template image, to allow comparisons between subjects. As for DW-MRI 
data, normalized and corrected functional scans were sub-sampled to the 
anatomical labeled template of the human brain [12]. Regional time series 
were estimated for each individual by averaging the fMRI time series over all 
voxels in each of the 90 regions. To eliminate low frequency noise (e.g. slow 
scanner drifts) and higher frequency artifacts from cardiac and respiratory 
oscillations, time-series were digitally filtered with a finite impulse response 
(FIR) filter with zero-phase distortion (bandwidth 0.01 — 0.1 Hz) as in pT3| . 

A functional link between two time series xfi) and Xj(t) (normalized 
to zero mean and unit variance) was defined by means of the linear cross¬ 
correlation coefficient computed as r %3 = (xi(t)xj(t)) , where (•) denotes the 
temporal average. For the sake of simplicity, we only considered here cor¬ 
relations at zero lag. To determine the probability that correlation values 
are significantly higher than what is expected from independent time series, 
Tij( 0) values (denoted r t] ) were firstly variance-stabilized by applying the 
Fisher’s Z transform 

Z , 3=0 .51n(l±W) (1) 

Under the hypothesis of independence, Zi 3 has a normal distribution with 
expected value 0 and variance 1 /{dfij — 3), where df is the effective number 
of degrees of freedom [Tii nil m- If the time series consist of independent 
measurements, df\j simply equals the sample size, N. Nevertheless, autocor- 
related time series do not meet the assumption of independence required by 
the standard significance test, yielding a greater Type I error da nans]. 

1 http://www.fil.ion.ucl.ac.uk 
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In presence of auto-correlated time series df must be corrected by the fol¬ 
lowing approximation ~ ^ ^ where r xx (r) is the 

autocorrelation of signal x at lag r. 

2.2 Networks Normalisation 

From Eq. (JT|) our networks weights are in a non-normalised interval Z l3 G 
[a, b] C R. In order to apply the framework described in [TS], we normalise 
our networks weights into the unit interval I = [0,1] by means of a unique 
linear function: 


(1 — 2e)Zij + (2e — 1) • MIN(Zij) 
MAX(Zij) - MIN(Zij) +e 


( 2 ) 


where e in general is set to 0.01 in order to avoid merging and isolate vertices 
with weights at the boundaries of Z t3 G [a,6]. As proved in [TS], since the 
normalisation is done by a unique linear function this does not affect networks 
properties. 

2.3 fMRI Networks aggregation and embedding 

Among many ways to aggregate a group of networks here we employed a 
topological algebraic way to aggregate a group of networks. The networks 
group possess the same nodes but different edges values and can mathemat¬ 
ical be represented by a weighted graph G = (N,E). N is the set of nodes 
representing the brain ROI’s (N = 90 in this study) and E is the set of 
edges values (connections) between ROI’s, e.g. G E : G [0,1] in the 

proximity space or Vdjj G E : dij G [0, +oo] in the distance space. 

For mathematical notation simplicity, we denote a network with the same 
notation we use to the set of nodes N. That is, a set of n networks (e.g. group 
of subjects) is represented by N k with k G {1,2,3...., n). 

One possible way to aggregate a group of n networks is simply by aver¬ 
aging the homologous edges values. Obtaining in this way a group represen¬ 
tative network, N*. 



m 



( 3 ) 


\k] 

where ejj is the edge e t j from network N k . 
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Another way to aggregate networks as explained in Simas et. al. [19], is 
by considering all networks as a multilayer network (often called multiplex ), 
which can be represented as fourth-order tensor |19j . This tensor can be 
represented as a extended matrix |21j . The work of Simas and Rocha [18] , 
introduces a framework to aggregate networks in an algebraic way, relating it 
with fuzzy logic reasoning, and in [T9] this work was extended for multilayer 
networks. In order to work algebraically with networks we have to set an 
algebra (defined as a vector space equipped with a bilinear product). This 
algebra allows us to perform algebraic operations with networks in the same 
way we perform algebraic operations in other contexts with other algebras 
(such as adding and multiplying real numbers). In short a network can be 
represented by an adjacency matrix and a multilayer network by a tensor. 
Considering a set of tensors working under the algebra L = (/, ©, ©), where 
the weights (tensor entries) of the tensors in / C R (subset of extended 
real line) and © and © two binary operators we may represent a multilayer 
network with tensor T in this algebra. In Simas et al. [ID] we have seen for the 
particular case of multiplex networks where layers are connected with weights 
Wi } i,L k ,L 3 = 1 (hr the proximity space), the representative group network (e.g. 
functional) can be represented by N* in the distance space (see below and 
Eq. [b]), as: 


N* = TV, © N 2 © • • • © N k (4) 

and the respective embedding by the following equation: 

Nembedded = N* © N* 2 © • • • © N* r (5) 

where N* is defined in Eq. Q and r, is the convergence parameter [IB] HD] - 
Figure [I] summarises the metric embedding of a multiplex network described 
above. 

Embedding a network of networks or, in our specific case, a a multiplex 
fMRI network, allows us to determine which edges in the several layers con¬ 
tribute to the aggregation. We can therefore determine the subjects that 
contribute more/less or none to the aggregated network, and identify in each 
subject the sub-graphs for which they contribute more. 

For our particular case, we embed our networks using the Metric Closure 
[18] defined by the algebra L = ([0, +oo], min, +), where © = min and 
© = +. The metric closure or metric embedding of a given network into a 
metric space, is a generalisation of All Pairs Shortest Paths Problem (APSP) 
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Multiplex Brain Networks 


Algebra 

L = (TV, ©, ®) 



Figure 1: Schematic representation of the main steps for the described net¬ 
works aggregation and metric embedding (defined here for the algebra L). 


as shown in Ref. [18], e.g. Johnson Algorithm can be used to calculate the 
metric closure [20] . 

Note that to calculate the metric closure (or Johnson algorithm) of a 
network we have to translate our networks from a proximity space into a 
distance space. There are many possible mappings to map a similarity space 
into a distance space, see [18]. Applying Eq. (J6]) to all network weights, 
Wij G [0,1] (for more details see [IS]), we obtain the isomorphic distance 
network with weights d itJ e [0, +oo]. 

dij = — - 1 ( 6 ) 

w id 

Embedding networks or multilayer networks allows us: (a) to detected 
clusters of nodes in a high-dimensional topological spaces, and by projecting 
the algebraic high-dimensional embedding into 3D, it allows to perform ex- 
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ploratory networks analysis and; (b) to preserve the multilayer sub-structures 
across layers/subjects, better than other aggregations methods, as we show 
for the specific case of ” simple” averaging (equation [3]) . 

Next we compared both methods of aggregation, ’’simple” averaging (eq. 
|3| and algebraic aggregation (eq. [5] according to [TO]) of our fMRI networks 
respectively, using our proposed method of embedding and comparison net¬ 
works. 

2.4 Multimodal Networks comparison 

In general networks have been compared using statistical measures of lo¬ 
cal and global properties of networks, such as: clustering coefficient, small¬ 
worldness, degree distributions, etc. We can fold in the literature some exam¬ 
ples of such techniques to compare multiple networks Our approach 

in this work is different. After embedding networks into the same metric 
space defined by the applied algebra, in our case L = ([0, +oo], min, +), we 
are able to compared them topologically. However, since networks generally 
come from different modalities (e.g. fMRI and DTI) it requires a previous 
step. We need to normalise the embedded edge weights distributions from 
the different modalities to the same average and variance to remove scale 
factors. One possible way to normalise both distributions, if we assume nor¬ 
mality, is by calculating the z-score of the edge weights distributions (zero 
average and standard deviation set to the unit). 

The embedded networks represent a hyper-grid in a multi-dimensional 
space with dimension equal or below to the number of nodes. In order to 
simplify and have some visual insight we can downgrade linearly this multidi¬ 
mensional grid into a 3D grid. This can be achieved applying to the embed¬ 
ded networks any technique for dimensionality reduction such as linear/non¬ 
linear Multi-Dimensional Scaling (MDS). MDS procedures refer to a set of 
related ordination techniques used in information visualisation, in particular 
to display the information contained in a distance matrix [23]- These tech¬ 
niques guarantee, with a given distortion, that the relative distance between 
nodes is preserved in both multi-dimensional and low-dimensional reduction 
space. Plotting this low-dimensional grid (e.g. in 3D) we can use any statis¬ 
tical technique to fit a continuous surface into the data (see below Figs. (|3j 
0 )- Its is natural to think that the difference between two surfaces obtained 
from different networks will emphasise topologically differences between the 
two connectivities. In this work we performed this operation in the multi- 


dimensional space by subtracting homologous embedded edges weights and 
take the absolute value of both embedded hyper-grids. This is, we subtract 
homologous embedded edges pairwise according to the formula: 
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M is the difference grid in the multi-dimensional space. Because the M-grid 
represents the difference between the two grids from different modalities (see 
above), the relative distance between nodes in M (given by Eq. ([?])) should 
be concentrated at the origin if they are topological similar, otherwise widely 
distributed in the multi-dimensional space. Nodes at a distance from the 
origin of k x a are statistical different of k standard deviations. Moreover, 
since we z-scored both embedded edge distributions this give us some degree 
of statistical significance when we compare both networks. All nodes that 
relay outside of a hyper-sphere with centre at the origin with radius R = kxa, 
are statistical different. Here we had set er = 1 for both distributions (z-score 
variables are estimated from the distributions of the embedded weights). 

Figure [2] illustrates this process. After applying to both fMRI* and DTI 
networks the same algebra and the metric embedding described above, both 
networks rely on the same metric space, therefore comparable. Topological 
differences can be visually seen in a linearly downgraded to 3D dimensions 
using a multi-dimensional scaling technique, which preserves the relative dis¬ 
tance between points in the grid (nodes or brain areas). 


3 Results 

In figure [3] we illustrate the results of different aggregation procedures on 
the ensemble of fMRI networks. Compared with a fMRI connectivity matrix 
from a single subject (Fig. |3](A)), one can notice the difference of a single av¬ 
eraging across subjects (Fig. [3](B) ) and our proposed algebraic topologically 
aggregated connectivity network (Fig. [3](C)). It is clear that the averaging 
procedure tends to blur connectivity values between nodes. In contrast, the 
topologically algebraic aggregation can preserve components that are com¬ 
mon across subjects. As other multilinear algebra or tensor-based analysis, 
our approach provide a natural mathematical framework for studying con¬ 
nectivity data with multidimensional structure. For illustrative purposes, we 
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Figure 2: Topological algebraic networks comparison. Connectivity from 
different modalities (here fMRI and DTI) are firstly embedded (black dot 
points on the manifolds indicate the brain nodes) and then compared in a 
low-dimensional space. Black points outside the sphere correspond to nodes 
with a topological difference (at a given threshold) in the two modalities. 


also show the DTI connectivity matrix in figure |3](D)). It worths noticing 
the similarity of the anatomical connectivity structure with the aggregated 
(multiplex) connectivity obtained in figure [3](C). Moreover, since each layer 
encodes the functional network for a given subject, each subject contributes 
to the tensor aggregation/embedding with some or none connections (edges), 
as depicted in metric closure, Eq. (§■ If a layer do not contribute for the 
aggregation/embedding, we may consider this layer (subject network) as an 
outlier. Moreover, we are also able to identify the specific sub-network con- 
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Figure 3: (A) fMRI single subject network (B) Average aggregated fMRI 
network (C) fMRI Algebraic Topologically aggregated (multiplex) network 
(D) DTI network. 


tribution (edges) of a given layer to the aggregation/embedding. 

Low-dimensional embeddings of different aggregated networks are illus¬ 
trated in figure [4j High-dimensional data, such as the information contained 
in the distance matrix obtained for the different networks, can be difficult 
to interpret. Here, multidimensional scaling (MDS) was used for visualising 
the level of similarity of individual nodes of each -aggregated- network. The 
MDS algorithm aims to place each node in a low dimensional space such that 
the between-nodes distances are preserved as well as possible. This repre¬ 
sentation into a low-dimensional space enables an exploratory analysis and 
makes data analysis algorithms more efficient. Indeed, from the different 
plots of figure [4] one can identify brain areas that are topologically close in 
the aggregated network as those points that are close on the 3D grid. This 
is clearly illustrated by the MDS representation of the multiplex functional 
network (Figure [4] (C)). Nodes from the occipital regions form a compact 
group of nodes topologically close (with similar connectivity structure), as 
revealed by the blue points depicted on Figure [5] (C). We also notice that a 
compact group of nodes is formed by regions of the temporal lobe, putamen 
and insula, which are indicated by the red circle. Similarly, the anatomical 
network in Figure |4](D) clearly displays a natural organisation, i.e. nodes of 
the two hemispheres lie on both sides of the dotted black line. Further, nodes 
from occipital regions in the anatomical network, indicated by the blue cir¬ 
cles (including calcarina, cuneus, precuneus, ...), are distantly located from 
the group of frontal brain areas indicated by the red marks. 

Finally, figure [5] displays the difference grid M in a low-dimensional space. 
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Figure 4: Multi-Dimentional Scaling (MDS) of the embedded networks (a) 
fMRI single subject (b) fMRI average embedded network (c) fMRI Algebraic 
Topological aggregation (multiplex) embedded network (d) DTI embedded 
network. Black dots indicate the embedded nodes. In plots (C-D), blue and 
red points indicate the groups of brain areas discussed in the text. 

As defined in Eq. ( [T]) , M corresponds to the relative distance between nodes 
in networks from different modalities. Differences between brain areas are 
represented as points widely distributed in the low-dimensional space. Those 
nodes from different modalities (fMRI and DTI) that share an identical topo¬ 
logical structure are located at the origin. The larger the difference in the 
connectivity structure, the larger the distance from the origin. By setting 
a threshold k x a, one can identify brain areas with similar connectivity as 
those points that lie inside of the hyper-sphere of radius k x a with centre 
at the origin. 

The number of brain regions (ROIs) with similar anatomical and func- 
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Figure 5: Comparisons between DTI and all other embedded fMRI networks. 

(A) 3D projections from Eq. ([7]). Only points outside the sphere are plotted. 

(B) Number of ROI’s inside the sphere of radius of k x a. Results from a 
single subject, average connectivity and multiplex networks are represented 
by the red, blue and green points and curves, respectively. We consider the 
regions statistically different for k > 1 and statistically equal for k < 1. 
This shows that the multiplex algebraic aggregation (green) is more similar 
algebraically to DTI then average aggregation (blue) and single subject fMRI 
network (red). 


tional connectivity are given in figure [5] (B) as a function of the threshold 
k x cr. Curves correspond to the number of regions inside a hyper-sphere 
of various radius. We notice that the number of regions differ as a function 
of the aggregated network’s type. It is worthy to mention that the differ¬ 
ences above /c-standard deviations are the important ones, since is above this 
threshold that the ROI’s or nodes become statistical different when compare 
networks. In our example, the fluctuations below one standard deviations 
may give us some trend but all nodes in the networks are statistical equal 
for all types of aggregation. For our specific case, as an example, the brain 
areas located outside the hyper-sphere of radius k — 1.2 for the two types of 
aggregation, are listed in the table 1. 


4 Discussion 

The recent prevalence of applications involving multidimensional and multi¬ 
modal brain data has increased the demand for technical developments in the 


13 





Averaged aggregated network 


Calcarine L 

Calcarine R 

Cuneus L 

Cuneus R 

Occipital Inf 

Lingual L 

Lingual R 
Temporal Sup L 
Temporal Sup R 
Insula R 

Occipital Sup L 
Occipital Sup R 
Occipital Mid L 
Occipital Mid R 

Multiplex aggregated network 

Cingulum Post R 
Amygdala R 
Postcentral R 

Occipital Mid R 
Occipital Inf L 
Occipital Inf R 

Thalamus L 

Heschl R 
Temporal Sup R 


Table 1: ROIs with connectivity differences from DTI at 1.2 standard devi¬ 
ation 

analysis of such complex data. Indeed, the discrepancy between structural 
and functional brain connectivity is a current challenge for understanding 
general brain functioning. In this paper, we presented a method for charac¬ 
terising the correspondences between functional and anatomical connectivity. 
To summarise, the main steps of our method are: 

1. Metric network embedding: This procedure embed a group of con¬ 
nectivity graphs in a common space allowing straightforward compar¬ 
isons. In contrast with simple averaging of connectivity matrices, the 
topologically algebraic aggregation can preserve components that are 
common across different subjects or different neuroimaging modalities. 
This tensor-based aggregation allows enhancing the common underly¬ 
ing structures providing a natural mathematical framework for study¬ 
ing connectivity data with multidimensional structure. 

2. Multimodal Networks comparison: the differences between the embed¬ 
ded networks are calculated and represented in low-dimensional space. 
Multi-Dimensional Scaling simply enables to display the information 
contained in the resulting distance matrix allowing thus an exploratory 
analysis of the data. 

3. Detection of nodes (ROI’s) with different connectivities: from points 
widely distributed in the low-dimensional space one can detect brain 
nodes that share a similar topological structure as those points are 
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located close to the the origin. One can identify brain areas with the 
largest difference between anatomical and functional connectivity as 
those points located outside an imaginary hyper-sphere of a radius 
given by a threshold (Table 1) 

Our findings suggest that embedding a brain network on a metric space 
may reveal regions that are members of large areas or subsystems rather than 
regions with a specific role in information processing. This is clearly illus¬ 
trated for the anatomical network in figure [5] D, where frontal and occipital 
brain areas of both hemispheres are situated at distantly and located points 
of the space. Contrary to a classical averaging of connectivity matrices, the 
embedding of the multiplex functional network reveals brain areas that play 
a role in large brain system such as the occipital regions, known to be active 
when the subject is at wakeful rest. 

Although experimental evidence suggests that functionally linked brain 
regions have an underlying structural core, this relationship does not exhibit 
a simple one-to-one mapping H- These correspondences have also been 
investigated in specific subsystems, must of them focused on the default 
mode network (DMN), which is a group of brain regions that preferentially 
activate when individuals engage in internal tasks, i.e. when the subject is 
not focused on the outside world but the brain is at wakeful rest. Several 
studies report that the DMN exhibits a high overlap in its structural and 
functional connectivity [251E]. Nevertheless, strong discrepancies have been 
reported and strong functional links can be found between regions without 
direct structural linkages [9]. 

At a group level, one of the reasons for this discrepancy between struc¬ 
tural and functional connectivity has been suggested to be the functional 
variability across subjects [26, 9[[25j. Indeed, clinical studies have provided 
evidence for a large heterogeneity of the functional connectivity, particularly 
in groups of patients with brain disorders such as neuropsychiatric disorders, 
which strongly alters the structural-functional relationships [25]. Analytical 
tools are therefore required to account for this variability in order to enhance 
the common underlying network structure. 

Results suggest that averaged aggregation captures the general differ¬ 
ences in regions that play a role in visual, auditory and body self-awareness 
processes, but fails to identify in detail other specific areas across the sub¬ 
jects/groups. In table 1 we observed that the average aggregation essentially 
captures part of visual (calcarine, cuneus, lingual, occipital), auditory (supe- 
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rior temporal gyrus), and insula regions that are associated to visual process 
and body self-awareness. We expected in average that these regions to be 
highly activated during the scans. This is alined with the fact that rs-fMRI 
was acquired with closed eyes and the subjects have some auditory and body 
awareness. 

From the multiplex aggregation (or algebraic aggregation) shown in table 
1, we observed that besides capturing the well known visual (occipital areas), 
primary sensory cortex (postcentral) and auditory regions (Heschl gyrus, su¬ 
perior temporal, thalamus), this approach also captures some other network 
sub-structures involved in touch activation (postcentral gyrus, thalamus) and 
emotional state activations (amygdala, thalamus, posterior cingulate). This 
alines with onr claim that algebraic aggregation preserves better the multi¬ 
layer sub-structures across a group of subjects (multilayers) accounting for 
as much of the variability in the data as possible. 

Although we cannot definitively a one-to-one mapping of the structural 
and functional connectivity, we think that our method could provide new 
insights on the organisation of brain networks during diverse cognitive or 
pathological states. We therefore hope that our approach will foster more 
principled and successful analysis of multimodal brain connectivity datasets. 

For all the methods described in this article we provide the corresponding 
MATLAB software code. Data and code are freely available at the website 
https://sites.google.com/site/fr2eborn/download 
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